Optimal Decompositions of Barely Separable States 
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Abstract 

Two families of bipartite mixed quantum states are studied for which it is 
proved that the number of members in the optimal-decomposition ensemble 
— the ensemble realizing the entanglement of formation — is greater than 
' the rank of the mixed state. We find examples for which the number of states 

in this optimal ensemble can be larger than the rank by an arbitrarily large 
factor. In one case the proof relies on the fact that the partial transpose of 
d ' the mixed state has zero eigenvalues; in the other case the result arises from 

the properties of product bases that are completable only by embedding in a 
larger Hilbert space. 
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I. INTRODUCTION 



The work of recent years has given us an extensive understanding of the entanglement 
of pure bipartite quantum states. While there are still many open questions about the 
entanglement of finite collections of quantum states [0], a rather complete understanding 
of 'asymptotic' entanglement, that of a large number of copies of a quantum state, has 
emerged: a pure state is unentangled if and only if the state can be written in a product form 
I \E') = I "^b)- The single good quantitative measure of entanglement is E = S{pa) = 
S{TrB \ ^)(^' I), where S is the von Neumann entropy. And, a collection of bipartite pure 
states with total entanglement E can be reversibly interconverted into any other collection 
of pure states with entanglement E by purely local operations. 

However, despite much recent effort, we cannot claim to have such a complete under- 
standing of quantum entanglement for bipartite mixed states. Its characterization has a 
remarkably greater complexity and richness than the pure-state case: There is not, except 
in very simple cases, an unambiguous way to say if a mixed state is entangled or not. 
There is no single good quantitative measure of mixed-state entanglement. And, it seems 
that entanglement is irreversibly lost when one attempts to convert it from one mixed-state 
embodiment to another. 

Much of this difficulty can be traced to the basic fact that there is no single way of 
viewing a mixed quantum state as an ensemble of pure states. In fact, we know that there 
are infinitely many such representations, and we have previously noted that in general these 
pure-state ensembles exhibit entirely different entanglement properties. For example, the 
completely mixed state of two qubits is equally well described as an ensemble of product 
basis states (no entanglement) or as an ensemble of the four Bell states (all maximally 
entangled). 

It seems that no single measure of entanglement for mixed states is correct, but many 
different ones are useful depending on the situation. The ensemble decomposition of a 
mixed state with the maximal entanglement is useful in situations where the two parties 



2 



holding the state can be given aid by a third party to extract a pure state with the greatest 
entanglement; we have termed the average pure-state entanglement of this ensemble the 
"entanglement of assistance" 0. Another, operational characterization of entanglement is 
the "distillable entanglement" D the average number of maximally entangled singlet 
pairs that can be extracted from many copies of the mixed state by local operations and 
classical communication. Yet another way of quantifying entanglement related to D has 
been proposed in which the minimum distance from the set of separable mixed states is 
taken as the measure of entanglement. 

Finally, the entanglement measure on which we focus in this paper is the "entanglement 
of formation" 0,0. It is the average pure-state entanglement of the ensemble which has 
minimal entanglement that describes the mixed state. Thus this is dual, in an operational 
sense, to the entanglement of assistance. It plays several other roles as well: it is converse, in 
some sense, to the distillable entanglement, in that it gives the number of EPR singlet pairs 
needed to create the mixed state by local operations. It, like the measure of entanglement 
in Ref. |]^, provides bounds on the distillable entanglement, and thus on other quantities 
such as the quantum capacity of noisy channels that are of great current interest in quantum 
information theory. 

Thus, we believe that a complete characterization of the mathematical properties of the 
entanglement of formation should be valuable in the continued development of quantum 
information theory. In this paper we give new results on one particular feature of the entan- 
glement of formation, the least number of states needed to make up a minimal-entanglement 
ensemble of a mixed state. (In such optimal decompositions of mixed quantum states 
have already been considered, but with respect to a function related to, but different from, 
the entanglement of formation.) Determining such optimal decompositions gives informa- 
tion about the minimal-complexity procedures for creating a mixed state from a supply of 
EPR singlets. But besides the operational significance of our results, we believe that the 
characterizations we provide are of greater significance on account of the light they shed 
on the complexity and richness of the mathematical structure of this important concept in 
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quantum information theory. 

II. ENTANGLEMENT BASICS 

Let p be a density matrix on the bipartite Hilbert space Ha^'Hb and let Sp = {pi, \ ipi)} 
with Pi > be an ensemble into which p can be decomposed: 

P = I]P»l^»)(V'i|- (1) 

i 

The entanglement of formation of p is defined as [§] 

Eip) = mmY^pMi^^){^^\). (2) 
i 

where 

E{\ I) = ^(Tr^l I) = ^(Tr^l ^)(^ |), (3) 

where S{.) is the von Neumann entropy: 

5(p) = -Trplogp. (4) 

The minimization in Eq. @) makes an analytical computation of the entanglement of 
formation of mixed states a nontrivial task. Only in a bipartite Hilbert space 2® 2 has the 
problem of determining the entanglement of formation of any density matrix been completely 
solved, in the work of Wootters [|10|- Uhlmann p| has shown that every bipartite density 
matrix p admits an optimal decomposition, that is, the one that achieves the entanglement 
of formation E of p (Eq. (^), with at least TZ{p) and at most TZ{p)'^ different pure states, 
where 7?.(p) is the rank of the density matrix. We call the number of different pure states in 
an ensemble that forms a decomposition of a density matrix p the cardinality of the ensemble. 
We say that the optimal ensemble cardinality of a separable state p, which we denote by 
Cs{p), is k if at least k different pure states are required for a separable decomposition of p. 
Since the number of states must at least be sufficient to span the support of p. 
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Ce{p)>n{p), (5) 

directly giving Uhlmann's lower bound. However, it has been a open question whether there 
are states for which the optimal decomposition has more than 7^(p) different pure states. 
Note that the states in the decomposition of a density matrix p are always in the range of 
p. This means that if Cs{p) > 'R-ip) the states in the optimal decomposition will be linearly 
dependent. 

In this paper we present two sets of examples of separable bipartite states p for which we 
prove that the cardinality of the optimal decomposition of p exceeds TZ{p). These are the 
first examples of such states. Both types of examples can be found in principle in arbitrary 
high dimensions. 

It is useful to classify states according to their behavior under partial transposition. Let 
p'^^ = (1^ T)p where T is transposition of a matrix in a chosen basis, p is positive under 
partial transposition (PPT) if p"^^ is a density matrix, i.e., it has no negative eigenvalues. If p 
is negative under partial transposition (NPT) then p"^^ has at least one negative eigenvalue. 
It is known that for 2 (8> 2 and 2 (g) 3 systems, PPT is a necessary and sufficient condition for 
separability [|ll|]. For a bipartite state in a Hilbert space of arbitrary dimension, PPT is a 
necessary condition for separability 



III. SEPARABLE STATES AT THE BOUNDARY 

In this section we show that if a separable state and its partial transpose have unequal 
ranks, then one of them must have its optimal ensemble cardinality greater than its rank. 
From this result we prove that partial transposes of full-rank separable states that lie on 
the boundary of the set of PPT states have optimal ensemble cardinality greater than their 
ranks. Finally we give examples of such states for any n®n system. 

We start with the following straightforward observation: 

Lemma 1 Let p he a separable state on Ha ® 'Hb- Let p^^ = (1^ ® T)p where T is 
transposition in a chosen basis. Then 
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(6) 



Proof We prove that Cs{p^^) < Ce{p) and Ce{p^^) > Cs{p). Since p is separable, its 
optimal decomposition involves only product states: 

p= Pi\i^i){i^i\<^\4>i){4>i\- (7) 

i=l 

Then it follows that 

Ce{p) 

P^^= E P^\i^^){i^^\®\^:){^:l (8) 
i=l 

and thus Cslp"^^) is at most Cs{p). By performing the partial transpose again on p^^ we 
can prove the inequality in the other direction. □ 

We have seen that the optimal ensemble cardinality is invariant under partial transposi- 
tion. The rank of a density matrix p is not necessarily invariant under partial transposition. 
We can draw the following conclusion: 

Theorem 1 Let p be a separable state on Ha ^'Hb- If 

n{p)^n{p^-) , (9) 

then either p has the property that Cs^p) > TZ^p) or p^^ has the property that Cs^p^^) > 

7^(pTB). 

Proof This follows directly from Lemma |1| and Eq. (^. □ 

Where do we find separable states p with the property Eq. (|^)? For this we look at 
full-rank separable states that lie on the boundary of the set of PPT states. The following 
lemma, illustrated by Fig. [^, looks into this: 

Lemma 2 Let p be a separable state on T-Ca ^'Hb with full rank, TZ{p) = dimTi^ ®Ti.B- If 
p lies on the boundary of the set of PPT states, then 

n{p)>n{p^-). (10) 
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Proof The set of PPT states iSppT = {p \ p^^ > 0} is a closed convex set. For the 
separable states p on the boundary of this set, the state p'^^ has at least one eigenvalue 
which is zero, as p is arbitary close to entangled states pE for which p^^ has at least one 
negative eigenvalue. Thus p^^ does not have full rank, and for full-rank states p this implies 
7^(p) > 7^(p^B). □ 

One can remark the following: In 2 2 and 2 ® 3 all entangled density matrices have 
the property that p'^^ ^ |Tl||. Therefore any separable density matrix p that is on the 
boundary of the set of separable states and has full rank, fulfills the conditions of Lemma 
H. With Theorem |l] it follows that the partial transposes of these density matrices have the 
desired property, i.e., Cs^p^^) > lZ{p^^). 

A final comment about the results of this section: Eqs. and (^) show that for 
separable p, p = p^^ if all the states 0j are real in some local basis, so obviously the ranks of 
p and p^^ are equal. Thus, any state p satisfying Lemma must have complex state vectors 
in any separable decomposition in any local basis. Note that even real density matrices p 
sometimes have optimal decompositions which require complex vectors [|I^]. Readers may 
find it surprising that even for complex vectors, there exist sets | ipi) \ (pi) for which the set 
spans a space of a different dimension; but this is exactly the consequence of Lemma 



A. Examples 



The generalized Werner state [14| in n ® n is defined as 



(11) 



where | \1'^) = J27=o I maximally entangled state in n ® n and < / < 1. Let 

p(/) = Pw(/)'^^on7^„®7^„. (12) 



It has been shown by Horodecki and Horodecki that the state Pw(/) is separable for 
< f < 1/n. On the other hand, for l/ra < / < 1, p{f) is not positive semidefinite. 



Therefore the state pyf{l/n) lies at the PPT-NPT boundary. Upon inspection of Eq. ([TlD 
we see that the rank of p^{l/n) is full, TZ{p^{lln)) = v?. Thus we can use Theorem ^and 
Lemma ^ 

It is easy to show by direct calculation that p{l/n) has exactly n{n — l)/2 zero eigenvalues 
and hence its rank is n(n + 1)/2. Therefore by Theorem |I| and Lemma 0, we find that p(l/n) 
has an optimal ensemble cardinality of at least whereas the TZ{p{l/n)) = n{n + l)/2. For 
this state, the ratio ^^/^/-i^/"?? can be as large as 2 when n tends to oo. This ratio could 
be made higher if there exists a p in n ® n for which p"""^ can have more than n{n — l)/2 
eigenvalues; we have no indication that this is possible. But for the state p{l/n)®^, k tensor 
copies of the above state, the ratio of to TZ can be made arbitrarily large (going to 2^ for 
n oo). 



IV. BARELY COMPLETABLE SETS OF PRODUCT STATES 



In |]T5|,|T6[ the notions of an unextendible product basis and an uncompletable product 
basis were introduced. A product basis is a set of k separable orthogonal pure states in 
n®m. Considering the case k < nm, the basis is unextendible if there are no additional pure 
product states orthogonal to all the members of the basis; it is uncompletable (in n®m) if the 



number of such additional states is less than nm—k. In [|15| , pi3|l it is shown that the completely 
mixed state p on the Hibert space complementary to the space spanned by the unextendible 
product basis is entangled, but has the property that p^^ is positive semidefinite. In the 
following, we will need the notion of a local extension of a bipartite Hilbert space H, = 
Ti-A ® 'Hb- a local extension of 7-^ is a Hilbert space H' = {Ha © H'j^) ® {Hb © H'b)- 



In [jT5[ and |T6| the notions of an unextendible product basis and an uncompletable 



product basis were introduced [|T^. It was shown how to construct, from an unextendible 
product basis, a bipartite entangled state p for which p'^^ jg positive semidefinite. We will 
need the notion of a local extension of a bipartite Hilbert space Ti = Ha ® Hb'- A local 
extension of 7i is a Hilbert space H' = {Ha © H'j^) ® (T^s © H'b)- 
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Theorem 2 Let { \ ai) ® \ A)}'^^';^ be a partial product basis S in Ti = Ha ® Ti-B- If S is 
uncompletable in 7i, but S is completable in some local extension 7i' of 7i, then ps defined 



as 



1 



= dimn-\S\ I ^ ~ ^ ' ^""^ ' ® ' ' I ' ^^^^ 



=1 



has the property that 

Ceips) > n{ps). (14) 



Proof As the set of states S is completable in a local extension of 7i, the state ps is 



separable, by Lemma 2 of [15]. The idea of this Lemma 2 is that the completion of the set 
{| aj) ® I l3i)}\Ji in H' give rise to a separable state 



1 



Ps = 1^1' - E I I ® I A) (A I j ■ (15) 

But Ps is obtained from p'g by local projections on to Ha and Hb and therefore p5 is 
separable as well. However, since ps is uncompletable in Ti, ps cannot be represented as an 
ensemble of orthogonal product states of cardinality lZ{ps)- Thus any optimal decomposition 
of Ps must use non-orthogonal product states. The von Neumann entropy of ps is equal to 
S{ps) = ^ogTZ^ps) as Ps is the identity on a space of dimension lZ{ps)- In order to achieve 
this entropy, the optimal decomposition of ps has to use more than 7l{ps) product states, 
or 

CeiPs) > nPs), (16) 

because any density matrix p which is a mixture of only n non-orthogonal states has entropy 
strictly less than logn bits. □ 
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A. An Example 



In [|T3| an example was given of a set of orthogonal product states on 3 (8> 4 that is not 
completable in 3 4, but is completable in 3 5. We reproduce the states here: Consider 
the states Vi ® Wi, i = 0, . . . , 4 with Vi defined as 

= A^(cos— ,sin— i = 0, ...,4, (17) 

with h = + and N = 2/^/5~+^. The states Wj are defined as 

Wj = A^'(-y/cos(7r/5) cos(2j7r/5), y^cos(7r/5) sin(2j7r/5), 

y^cos(27r/5) cos(4j7r/5), y^cos(27r/5) sin(4j7r/5)), (18) 

with normalization A^' = y27\7^. Note that wjwj+i = (addition mod 5). One can 
show that this set, albeit extendible on 3 4, is not completable: One can at most add 
three vectors like vo ® (wq, w^i, ^^4)"'", ® (^2, "U^s, 1(74)-'- and {vo,V3)-^ (t^i, 1^2, 'U^4)-'-. The 
completion of this set in 3 5 is particularly simple, being given by the following ten states: 

{vi, Vi)-^ ® Xo, (w^ G span(x4, Xi)), 
(^^0, "^^2)"^ ®xi, vi® {wf- G span(xo, X2)), 

{vi,V3)-^ ® X2, 1/2 (S) (w^ G span(xi, X3)), (19) 

("52, V4)-^ (g) fs, Vs (g) {id^ G span(f2, X4)), 

(■i^o, ■ws)"^ ® X4, V4 (g) (wf G span(f3, xq)). 
The state on 3 ® 4 has rank seven, but the separable decomposition consists of ten 
non-orthogonal states obtained by projecting the orthogonal states of the completion, Eq. 
(p!9|), back into the 3 (8> 4 Hilbert space. It is not known whether there exists a separable 
decomposition with more than seven but with fewer than ten states. 

V. DISCUSSION 

The results presented here on the minimum cardinality of optimal ensembles raises a 
large number of tantalizing questions; we would like to know this cardinality for all possible 
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mixed states. So far, our rigorous results apply only to separable states. There is some 
empirical evidence that for entangled mixed states as well (in fact, for states arising in the 
theory of unextendible product bases), the optimal ensembles can have a cardinality greater 



than the rank [jT8[. But we have found no techniques for proving any results for inseparable 
mixed states. We would also like to know whether there are cases for which the Uhlmann 
upper bound of 7?.(p)^ is attained. The states shown above are still far from this; for the 



states at the end of Sec. |1II A| with n = 2 and any k, Ce = TZ"^ . Finally we note that all the 
rigorous results we have pertain to cases where Cs, while greater than TZ, never exceeds the 
Hilbert space dimension. Is there some reason that Cs can never exceed this dimension? In 
conclusion, we have shown two different families of unentangled mixed states for which the 
minimal number of pure states in an optimal minimal-entanglement ensemble is provably 
greater than the rank of the mixed state. In both cases the proofs are possible because the 
mixed state is marginally separable, in the first case because the partial transpose of the 
state has zero eigenvalues, and in the second because the state is defined as the complement 
of a barely completable product basis. 



ACKNOWLEDGMENTS 

We are grateful to Peter Shor for the construction of the barely uncompletable product- 
basis example used here. We thank Charles Bennett, Richard Jozsa, John Smolin, and 
Armin Uhlmann for helpful discussions. AVT would like to thank David D. Awschalom 
for his invaluable support without which it would have been impossible to work in this 
exciting field. This work has been supported in part by the Army Research Office under 
contract numbers DAAG55-98-C-0041 and DAAG55-98-1-0366. AVT and BMT would hke 
to thank IBM Research for logistical support during their visits to the IBM Thomas J. 
Watson Research Center. 



11 



REFERENCES 



[1] See, for example, D. Jonathan and M. B. Plenio, "Entanglement-assisted local manip- 
ulation of pure quantum states," |quant-ph /9905071 . 



[2] E. Schrodinger, "Probability relations between separated systems," Proc. Camb. Phil. 
Soc. 31, 446 (1935). 

[3] L. P. Hughston, R. Jozsa, and W. K. Wootters, "A complete classification of quantum 
ensembles having a given density matrix," Phys. Lett. A 183, 14 (1993). 

[4] O. Cohen, Phys. Rev. Lett 80, 2493 (1998); D. P. DiVincenzo, C. A. Fuchs, J. A. Smolin, 
A. Thapliyal, and A. Uhlmann, "Entanglement of assistance", "Proceedings of the 
First NASA International Conference on Quantum Computing and Quantum Commu- 
nications", 17-20th February 1998, Palm Springs, CA, Springer- Verlag Lecture Notes 
in Computer Science Vol. 1509, (Springer- Verlag, Heidelberg, Germany, 1999) [LANL 
preprint |quant-ph/980303^ . 



[5] C. H. Bennett, H. J. Bernstein, S. Popescu, B. Schumacher, Phys. Rev. A 53, 2046 
(1996) [LANL preprint |quant-ph/9511030|l . 



[6] C. H. Bennett, D. P. DiVincenzo, J. A. Smolin, W. K. Wootters, Phys. Rev. A 54, 3824 



(1996) [LANL preprint |quant-ph/9604024|| . 



[7] V. Vedral, M. B. Plenio, M. A. Rippin, P. L. Knight, Phys. Rev. Lett. 78, 2275 (1997) 
[LANL preprint [quant -ph/ 9702027|| . A more detailed account can be found in V. Vedral 



and M. Plenio, Phys. Rev. A 57, 1619 (1998) [LANL preprint [quant -ph/ 9707035 
[8] A. Uhlmann, Open Sys. & Information Dyn. 5, 209 (1998). 

[9] It has been pointed out that a more satisfactory definition of the entanglement of for- 
mation would be 

E(p)= min lJ2P^Ei\A){i^^\). (20) 

N,tp^N iV . 

12 



Eqs. pop and are equal if the entanglement of formation is additive, but this property 
has not been proved. The analysis of this paper concerns the original definition of 
entanglement of formation, Eq. (Q). 



[10] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998) [LANL preprint [quant -ph/9709029| 



[11] R. Horodecki, M. Horodecki, and P. Horodecki, Phys. Lett. A 223, 1 (1996). 
[12] A. Peres, Phys. Rev. Lett. 76, 1413 (1996), [LANL preprint |quant-ph/9604005 



[13] Eq.(30) of [^] gives a decomposition with eight members for any mixture of Bell states. 
However it can be easily seen that appropriate sets of four of these give optimal- 
cardinality ensembles with at least one complex vector. 

[14] M. Horodecki and P. Horodecki, "Reduction criterion of separability and limits for a 
class of protocols of entanglement distillation," [LANL preprint [quant -ph/97080"l5[ . 



[15] C. H. Bennett, D. P. DiVincenzo, T. Mor, P. W. Shor, J. A. Smolin and B. M. Ter- 
hal, "Unextendible Product Bases and Bound Entanglement," [LANL preprint 
quant -ph/ 9808030. 



[16] C. H. Bennett, D. P. DiVincenzo, T. Mor, P. W. Shor, J. A. Smolin and B. M. Terhal, 
"Unextendible Product Bases and Bound Entanglement II," in preparation. 

[17] A product basis is a set of k separable orthogonal pure states in n^m. Considering the 
case k < nm, the basis is unextendible if these are no additional pure product states 
orthogonal to all the members of the basis; it is uncompletable (in n ® m) if the number 
of such additional states is less than nm — k. 

[18] J. A. Smolin (private communication). 



13 



FIGURES 




FIG. 1. Illustration of the construction of Lemma |2[ The PPT, NPT, separable, and inseparable 
sets of density matrices are indicated. Bi is the boundary between separable and inseparable, B2 
between PPT and NPT. For 2 2 and 2 3, Bi = B2, but in general they do not coincide. The 
conditions of the Lemma are met when p has full rank, is separable, and is on the boundary B2 
so that its partial transpose p^^ has some zero eigenvalues, p must therefore be at a place where 
Bi and B2 coincide, such as the arc ab. p"^^ is separable; also, it cannot be on the boundary B2, 
because its partial transpose, which is p, has strictly positive eigenvalues. However, it is possible 
that could sit on the boundary Bi where Bi and B2 do not coincide. 
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